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Abstract 

We study three classes of continuous time Markov processes (inclusion process, exclusion 
process, independent walkers) and a family of interacting diffusions (Brownian energy process). 
For each model we define a boundary driven process which is obtained by placing the system in 
contact with proper reservoirs, working at different particle densities or different temperatures. 
We show that all the models are exactly solvable by duality, using a dual process with absorbing 
boundaries. The solution does also apply to the so-called thermalization limit in which particles 
or energy is instantaneously redistributed among sites. 

The results shows that duality is a versatile tool for analyzing stochastic models of transport, 
while the analysis in the literature has been so far limited to particular instances. Long-range 
correlations naturally emerge as a result of the interaction of dual particles at the microscopic 
level and the explicit computations of covariances match, in the scaling limit, the predictions of 
the macroscopic fluctuation theory. 



1 Introduction 

Interacting particle systems are classical models to study non-equilibrium statistical mechanics. The 
standard setting is the one in which a system is placed in contact with reservoirs working at different 
parameters that create a stationary state characterized by a non-zero averaged current. The prototyp- 
ical example are the Symmetric Exclusion process with at most one particle per site connected to birth 
and death process at the boundaries [Lj |D] and the KMP process |KMP| connected to reservoirs which 
impose at the boundaries Boltzmann-Gibbs distribution with different temperatures. The Symmetric 
exclusion process is a model for transport of a discrete quantity, whereas the KMP process models 
transport of a continuous quantity. 

Problems that are very hard for classical Hamiltonian systems - for instance deriving Fourier law 
starting from the microscopic evolution - can be successfully approached using stochastic models. 
Furthermore stochastic models of transport have been used to prove new theorems in non-equilibrium 
statistical mechanics, such as the fluctuation theorem |GCllK] . to introduce new principles, such as the 
additivity principle jBDj . to construct new schemes, such as the macroscopic fluctuation theory that 
describes the density and current large deviations for diffusive systems [BDGJLll IDLSj , to test new 
algorithms, such as cloning algorithms to simulate rare events ( jGKLTj ) . Recently, the connection 
between deterministic Hamiltonian systems and stochastic models is emerging either by considering 
evolutions in which they are coupled |B0] or by considering slow/fast variables |DL) . 
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In the analysis of general stochastic processes a powerful tool, when available, is duality theory [51 Hi]. 
In the duality approach, the original process is studied via a simpler process, called dual process. This 
has been applied in different contexts, including interacting particles systems, interacting diffusions, 
queueing theory and mathematical population genetics (for a recent review on duality, which also 
include many references, see |JKj ). For stochastic models of transport, in the presence of boundary 
reservoirs, duality has been found originally for the boundary driven Symmetric Exclusion process 
with at most one particle per site |Spo| and for the KMP model |KMP| . Recently duality has been 
obtained also for other models, such as Symmetric Inclusion process and for interacting diffusion 
introduced to model heat conduction }GKR[ IGRV| . Typically the duality approach allows to study 
n-point correlation functions in the original process with L degrees of freedom using n interacting dual 
particles (with n << L). Thus the problem of studying correlations functions in a large (eventually 
infinite) systems is reduced to a finite dimensional problem. In all those cases the existence of a dual 
process has been related to the study of hidden symmetries and/or change of representation of an 
abstract operator [GKRVilSS] , 

In this paper we define classes of boundary driven interacting particle systems and a family of inter- 
acting diffusions with reservoirs. The particles systems interpolate from the Inclusion model with a 
Negative-Binomial stationary measure at equilibrium, to the Exclusion model with a Binomial equi- 
librium state, via Independent walkers with a Poisson stationary measure. The interacting diffusions 
are given by the so-called Brownian Energy process, having a Gamma distribution at equilibrium. 
For all the models we construct a dual process which allows to study correlation functions in the sta- 
tionary non-equilibrium state (i.e. different reservoir densities or temperatures) using the absorption 
probabilities of dual walkers at the boundaries. The results include, but are not limited to, some of 
the models mentioned above. The results hold true for any value of the reservoir parameters. This 
shows that duality is a versatile tool in the analysis of such systems. 

The KMP model is not included in the classes of model we consider in a first stage. However it can 
be recovered by considering an instantaneous thermalization limit of one element of the family of the 
Brownian Energy model, as already remarked in |GKRV| . We construct here the thermalization limit 
for all the models, obtaining both in the discrete and in the continuous case, several redistribution 
rules of particles/energy. A non-trivial stationary state is found in this case even considering only one 
site, since the reservoirs are not additive. 

Finally, we compute in this paper the explicit microscopic expressions for the averages and covarianccs 
for many parameters values of our models in non-equilibrium setting. From the analysis, we deduce 
that the macroscopic limit of the microscopic expressions are in full agreement with the predictions 
of the macroscopic fluctuation theory. 

The paper is organized as follows. In the next section we introduce the models and discuss the relation 
between them, including scaling limits. We describe the equilibrium states in section [21 both in the 
case of free boundary conditions and in the case of systems connected to reservoirs. Section [4] presents 
the duality result and describes the dual processes, together with the illustration of the use of duality 
in the non-equilibrium setting. As an example, the linear temperature profile are computed for all the 
models, and the stationary measure of boundary driven independent walkers is deduced from duality. 
Section[5]describes the instantaneous thermalization limit of our models, and their duals. The last two 
section deal with correlations in the stationary non equilibrium state (section ^ and the comparison 
to the predictions of the macroscopic fluctuation theory (section [7|). 

2 Models definition 

In this section we introduce our models. In the most complete setting, they are constituted by a 
bulk which is kept in a non equilibrium state by the contact with particles or energy reservoirs. In 
particular, we consider one-dimensional systems on a finite lattice {I, . . . whose boundaries (i.e. 
sites 1 and L) interact with the reservoirs. When needed, the reservoirs themselves will be represented 
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by two extra sites, namely sites and L + 1. 

Accordingly, the generators of the random processes associated with our models can be generically 
expressed as the sum of three terms 

C = Ca+Co+Cb, (2.1) 

where Cq represents the generator of the dynamics in the bulk, while Ca and Ch represent the gener- 
ators of the reservoirs. 

We will consider four models: three classes of interacting particle systems, characterized by the dif- 
ferent interactions between the particles, and one family of interacting diffusions introduced to model 
heat conduction pCRl iGKRVl IGRV] . The models are: 

1. the Symmetric Inclusion Process (SIP), with attractive interaction between neighboring particles; 

2. the Symmetric Exclusion Process (SEP), with repulsive interaction between neighboring parti- 
cles; 

3. the Independent Random Walkers (IRW), without interactions among particles; 

4. the Brownian Energy Process (BEP). 

In the first three cases the dynamic variable is a vector that specifies the number of particle on each 
site: 77 = (771, ...,r?L) £ O; here fl, the state space, depends on the model and will be defined ahead. 
In the case of the BEP the dynamic variable is a vector z representing the energies on each site of the 
lattice: z = (zi, . . . ,zl) G ^ = M+. 

2.1 Interacting particle systems 

The generators of the reservoirs for SIP, SEP and IRW have the following general form: 

jOafiv) = b{m)[f{v°^')-f{v)]+d{vi)[f{v'-')-.f{v)], (2.2) 
= b{r^L)[f{rj'^+'-^)-firj)]+dir^i^)[f{rj'^-'^+')~m]. (2.3) 

According to (|2.2p and (j2.3[) particles are injected into the system through the boundaries with rate 
b{n) with 7T, G No, and removed from the same sites with rate d{n). While b{n) is model-dependent, 
the annihilation rate is not, being in any case proportional to the number of particles at the boundary 
site. 

We introduce now our models by defining the actions of the generators C on the functions / : $7 — > R. 
Hereafter, we will use the symbol to denote the configuration obtained from r/ by moving a 

particle from site i to site i + 1, i.e. = (^r]i, . . . , 77^ — 1, r/^+i 4- 1, . . . , ?/l)- 

Inclusion walkers SIP(2fc). This is an inclusion process in which each site can accomodate an 
arbitrary number of particles, thus il — Nq . In the bulk each particle may jump to its left or right 
neighbouring site with rates proportional to the number of particles in the departure site and to 
the number of particles in the arrival site. In each boundary site particles are created with a rate 
proportional to 2k plus the number of particles sitting in that site; k G M_|_ labels the class of models. 
The generator is 

= a{2k + ,n) [fiv"-') fiv)] + im - firi)] 

L-1 

i=l 

+ 5{2k + tjl) [/(7;^+1-^) - fiv)] + PVL [/(»y^^^+^) - fiv)] • 
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The positive numbers a and 7 (resp. S and /3) tune the creation and annihilation rates of the left 
(resp. right) reservoirs. 



Exclusion v^ralkers SEP(2j). For this process the maximum occupation number at each site is 
2j e N, thus f2 = {0, 1, . . . , 2j}^. In the bulk particles jump independently to nearest neighbouring 
lattices sites at rate proportional to the number of particles in the departure site times the number 
of holes in the arrival site. The reservoirs inject particles in the systems with a rate proportional to 
the holes in the boundary sites. The generator is 

C'^^m = + + £r^/(,y) (2.5) 

= a{2j - rn) [/(r?"'i) - fin)] + 7^ [fiv'-") - f{v)] 

+ m ~ Vl) [/(r/^+^'^) - fiv)] + PVL [/(r/^'^+^) - fiv)] • 
The parameters a, 7, S, (3 have the same meaning as in the SIP(2fc). 

Independent random walkers IRW. In this case each particle jumps independently to nearest 
neighbouring lattices sites at rate 1, and each site can accomodate an arbitrary number of particles, 
thus ri = Nq . Jumps occur with the same probability to the right and to the left, while particles 
are created at rates a and 5 irrespective of the numer of particles at the boundaries. Therefore the 
system is described by the generator 

C^'^'^m = /:f^/(^)+4«>^/(^)+£^«^/(,y) (2.6) 
= «[/(^y"-')-/(r/)]+7r/i[/(r?^'")-/(r/)] 

L-l 

1=1 

+ 5 [/(77^+i'^) - fill)] + [iilL [/(r/^'^+^) - .f{^)] . 

The dynamics in the bulk can be further described by saying that if at site i there are r]i particles, one 
of the particle jumps at rate 77,; either to the left or to the right. As in the previous cases, parameters 
7 and P define the annihilation processes. 

Remark 2.1. The effect of the reservoirs is to impose the average number of particles on the left and 
on the right sides of the chains. With some misuse of language, but sticking to standard notations, we 
will call "densities" these averages and we will denote them pa (left reservoir) and pf, (right reservoir) . 
The vaules of pa and p^ are reported in the table below and computed in 5ecO 



System 


Pa 


Pb 


SIP 


■y — a 




SEP 




'^■1 0+S 


IRW 


a 
7 


b 





Table 1: Definition of pa and pjj. 



Remark 2.2. Note that the SIP process requires 7 > a and (3 > 5. This condition turns out to be 
necessary in order for the system to reach a stationary state (see also formula l|3.8p ). 
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Remark 2.3. It is interesting to remark that the exclusion (resp. inclusion) walkers with parameters 
(a, 7', (5, /?') converges to the independent walkers with parameters (a,7,(5, /3) in the limit j — ^ 00 (resp. 
k ^ 00) under the scaling 7' = 2j^, /?' = 2j/3 (resp. 7' = 2kj, /3' = 2kp). Indeed, in this limit the 

rSIP ^SEP r>TnW 

generators — converge to C '^^ , This remark can he put on rigorous grounds by using 

the Trotter-Kunz theorem (see Theorem 2.12 of Wj): see for istance IGKRVf for the proof in the case 
of SEP(2j). 



2.2 Interacting diffusions 

The last process we consider, the Brownian Energy Process, was originahy introduced in |GKRV] 
and recognized as dual of a system of inclusion random walkers. Here wc present the boudary driven 
version of the BEP. 



Brownian energy process BEP(2fc). This model describes symmetric energy exchange between 
nearest neighbouring sites. The dynamical variables (energies) are collected in the vector z = 
(zi, . . . , zl) G M.^ and the generator is 

^BEP ^ C^^^fir^) + C^^Pfi,^) + ci^Pfir^) (2.7) 
_d_ 

" dzi ^^dzlJ 2^^dzi 



L-l / Q o \ 2 



d , ( d d 



7^ - - 2fc(z, - z,+i) — - 



d d''\ I d 



Remark 2.4. The origin of the bulk dynamics, generated by 



dzi dzi+i J \dzi dzi+i 

can be explained as follows \GKR[ \ GKRVf . Consider m = 4fc G N velocity variables on each site i and 
call them Vi^a with a = 1, . . . ,m. Suppose that they evolve with the following generator 

which defines a process, called Brownian Momentum Process, introduced in \B0[ \ GE^ . Each term in 
£BMP j-gpj-gsg7jis a rotation in the plane {vi^onVi+i^/i) , therefore it conserves the total length vf ^ + 
'^1+1 B' '-^^ total kinetic energy. One can check that the BEP{2k) is the evolution process, induced 
by i2.9\) . of the total energies on each site 

m 

z, = Y,vl^. (2.10) 

a = l 

The generator of the BEP reservoirs C^^^ and C^^^ , that will be discussed in some details in Sec. 
impose an average energy AkTa on the left, and an average energy AkT}, on the right. The choice 
of their form is motivated as follows. Consider an Ornestein- Uhlenbeck process on each of the m 
velocities at site 1 of the Brownian Momentum process i2.9\) . namely 



^BMP ^ ^ 2T-V - vi,c.jr— ■ (2.11) 

a— 1 ' 



92 d 
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Since in the stationary state of this reservoir the {vi^a}a=i,...,m o-fe independent centered Gaussian 
with variance T then, using i2.10\) . the expectation of zi is E(zi) = X]™=i ^^(''^i a) ~ "^"^ ~ 4A:T . 

2.3 Scaling limit of the particle systems 

Besides duality, there is another relation eonnecting the bulk part of the BEP with generator Cq^^ 
(HiH]), and the bulk part of the SIP with generator C^'^ (third line in (EH)). The BEP can be 
indeed obtained from the SIP, through a suitable scaling limit, by a reinterpretation of this process 
as a model of energy transport, by supposing that each particle carries a quantum of energy e. In 
this interpretation, since Cq^^ conserves the number of particles, then it conserves the total energy. 
Consider the free boundary inclusion process r]{t) = (771 (i), . . . ,r]L{t)) generated by Cq^^ and let N 
be the total number of particles, i.e. = X^i^i Vi- Let e be a parameter of the order of 1/N, then one 
expects r]i to be of the order of as e — >^ (despite attractive interactions for any finite k there are 
no condensation phenomena in the SIP; one needs to rescale k with e to see particles coalescing into 
a single site; see jGRV2j ) . Then one may investigate the continuous dynamics generated in the limit 
as e — )■ on the variables Zi(t) ~ ei}i{t). It turns out that the limiting dynamics for z{t) is generated 
by CE'^r 

Proposition 2.5. Let 7]{t) = (?/i(0' • • • iVLit)) be the bulk inclusion process generated by Cq^^ with 
N particles. Let e = £/N for some fixed E > Q. Then the process z(t) = {ziit), . . . , ZL{t)) where 
Ziit) = erji{t) is, in the limit e — >■ 0, the bulk Brownian energy process generated by Cq^^ with total 
energy £ . 

Proof. Let F : Mf^ M., F = F{z) be a two times countinuously differentiable function, i.e. 
F e C'^iR^). Let = (z|, . . . , z|^) e R;^ be such that z'-fe G N^, then for any F as above, there 
exists / : ^- R, / = 7(77), such that 



Fizl,...,zi) :=/ 



(2.12) 



Let Cf) be the generator of the process z'^{t) induced by the SIP, then acts on F = F{z'^) as follows: 
[QF]iz^) = [C','^f] 



L-l 



e e 



1 ... 



/If 



+ 



2k- 



L-l 



+ 



2k 



/ 



I 



[F{zl,...,zt-e,zl+^+e,...,zl)^F{z% 
[F{zl,...,zt + e,zt+^-e,...,zl)-F{z')\ 



Suppose that converges to a finite limit — > z £ M;^ as e — > 0. Then, from the regularity 
assumptions on F , we have 



[^h+iF]{z') ~ Fizl . . . , zt_, -e,zt + €,..., zl) - F[zl, z|_i, . . . , z^) 

d d 



dzi 9z,+i 



F{z^) + o{e) 



(2.13) 
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while 



d d 



dz, dzi- 



F{z') + o{e), (2.14) 



and 

2 



Therefore we have 

d d \ , , ( d d 



F(z<^) + o(l). (2.15) 



Thus, for any F as above, lini£^o['Co-F](2:') — [Cq^^ F]{z). Moreover the total energy is clearly 
conserved in the limit and it is given by = z| = e iV ~ □ 

The same scaling analysis of the inclusion walkers can be performed on the bulk dynamics of inde- 
pendent random walkers. This yields a deterministic process as scaling limit, which is also dual to 
independent random walkers (cfr. jGKRVj ). 

Proposition 2.6. Letrj{t) = (771 (t), . . . , 7yi(t)) be the bulk process generated by C^^^^ with N particles. 
Let e = £/N for some fixed £ > 0. Then the process y(t) = {yi{t), . . . ,yL(t)) where yi{i) = erji{t) is, 
in the limit e 0, the deterministic energy process (DEP) with total energy yi(^) — ^ generated 

by 

^"^ ' d d 



dyt+i dyi 



Remark 2.7. One may wonder whether there exists a diffusion process arising as a limit of the 
Exclusion process. By performing an analogous scaling as above, the rates of the SEP take the form 
Nzi{2j — N Zi^i) that become negative in the limit as N ^ 00. Consistently the limit of the SEP 
generator is a second order differential operator that cannot be interpreted as the generator of a Markov 
process, since it is has a negative coefficient in front of the second order derivatives, i.e. 

L ( Q (?\^ f d d \ 
V -z,z,+i 7^ - 2j(z. - z,+i) — - . (2.16) 

Remark 2.8. The same scaling limit which transfoms the hulk dynamics of the SIP into the one of 
the BEP does not work with the reservoirs. Indeed, applying to C^^^ the scaling of Propositon \2.5l 
the resulting generator is 

(«-7)^i^ (2.17) 

which produces a deterministic behavior: zi{t) ~ zi(0)e'"~'''-'* . On the other hand it is simple to 
check that the thermal bath of the BEP can be obtained from a boundary driven SIP with a modified 
reservoir generated by 

= (^kq + (9 - [/(^°'') - fiv)] + qvi [fiv'-') - fiv)] (2.18) 
with the condition qe — )■ Ta as e ^ Q. 
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3 Stationary measures at equilibrium 



The models introduced in the previous section are Markov processes with discrete or continuous state 
spaces. The long term behaviour of the processes are described by their stationary measures. In 
general it is hard to determine such measures and, in fact, the invariant states of SIP, SEP and BEP 
in non-equilibrium conditions are not explicitly known. The problem of finding the explicit form of the 
invariant states is greatly simplyfied at equilibrium. The equilibrium condition for our systems can 
be obtained in two ways: either by suppressing the reservoirs (i.e. considering only the bulk dynamics 
Cq) or, retaining the reservoirs, by imposing equal densities or equal temperatures at the boundaries 
of the chain, i.e. p = pa = Pb or T = Ta = Ti,. 

In the first case there exists an infinite family of reversible measures labelled by a continuous parameter. 
In the second case (i.e. in the presence of the reservoir) at density p (resp. at temperature T) the 
boundary conditions select one reversible measure. 



3.1 Equilibrium product measures 

Reversible invariant probability measure P of the bulk dynamics generated by Cq can be obtained by 
imposing the detailed balance condition. When the state space is finite or countable, this condition 
is expressed by requiring that for any pair of configurations rj, ?/ G Q the probability P satisfies 

Lo{v,v')nv)^L„{r^',v)nv') (3.1) 

where L(){rj,r]') is the transition rate from the configuration rj to rj'^ i.e. Lo{ri,ri') ~ Cof{r]) with 
f{ri) — Sri^n'- When the state space 51 is continuos, a probability measure with density 'ip{x) is said 
to be reversible stationary measure if, for all functions / and g in the domain of the generator Cq, it 
holds 

f{x)Cog{x)^ix)dx ^ / Cofix)g{x)^{x)dx . (3.2) 



By imposing in the case of SIP, SEP, IRW and (j3.2|) in the case of BEP and requiring the 

factorization of the probability measure one obtain the reversible measures described in the following 
proposition, whose proof is left to the reader. 

Proposition 3.1. For the bulk processes with generator Cq defined in Sec. 2 we have 
Inclusion walkers SIP(2fc) 

The process with generator Cq^^ has a reversible stationary measure given by products of generalized 
Negative Binomial measures with parameters 2fc > and arbitrary < p < 1, i.e. 

Exclusion walkers SEP(2j) 

The process with generator Cq^^ has reversible stationary measure given by products of Binomial 
measures with parameters 2j G N and arbitrary < p < 1, i.e. 

Independent random walkers IRW 

The process with generator Cq^^ has reversible stationary measure given by products of Poisson 
distribution with arbitrary parameter X > i.e. 

m = t[^e-\ (3.5) 

i=i 
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Brownian energy process BEP(2fc) 

The process with generator Cq^^ has reversible measures given by product of Gamma distributions with 
parameters 2k > and arbitrary > 0, i.e. 

3.2 Equilibrium product measure with reservoirs 

We recall that, in the case of particle systems (see SecH]), the reservoirs are modeled by birth-death 
processes with creation rate b(n) and annihilation rate d{n), n the number of particles at the bound- 
ary. Each reservoir has, thus, its own reversible invariant probability measure, p(n), which satisfies 
the detailed balance condition b{n)p{n) = d{n + l)p{n + 1). This condition can be used to compute 
p{n). The average value of the random number n (that wc call density, irrespective to its value) is the 
quantity imposed by the reservoir to the system. 

The effects of the reservoirs, under the equilibrium conditions, arc described in the following proposi- 
tion, which can easily be proved with an explicit computation. 

Proposition 3.2. For the processes with generator C defined in Sec. 2 we have: 
Inclusion walkers SIP(2fc) 

The /e/t reservoir is modeled by the birth and death process with rates 

b{n)^a{2k + n), d{n) = -fn, neN. (3.7) 

The stationary state of this reservoir is given by a Negative Binomial measure with parameters 2k 
and p = ^. The reservoir density is p := (n) ~ 2k ~ 2k^:^^. The boundary driven process with 
generator C^^^ defined in {2.4^ , with parameters a, 7 and /?, 5 such that a/3 — -fS = (and thus 
Pa ^ Pb) admits the stationary product distribution: 

S 

Negative-Binom(2fc,p) with p := — = — for a < 7 and 5 < (3 (3-8) 

7 fi 

Exclusion walkers SEP (2j) 

The left reservoir is modeled by 

b{n)=a{2j ~n), d{n) = -fn, n G {0, 1, . . . , 2j}. (3.9) 

The stationary state of this reservoir is given by a Binomial measure with parameters 2j and p ~ ::^^- 
The reservoir density is p := (n) = 2jp = 2j;j^^. The boundary driven process with generator C^^^ 
defined in 12. 5\) . with parameters a, 7 and P, 6 such that a/3 — = (and thus Pa = Pb) admits the 
stationary product distribution: 

d S 

®ti Binom(2.?» with p = — - . (3.10) 

J + a p + 

Independent random walkers IRW 

The left reservoir has a constant birth rate 

b{n) = a, d{n) = jn, neN . (3.11) 

This reservoir imposes a Poisson measure with parameter X = ^. Therefore the density (i.e. mean 

number of particle) is p :~ {n) ~ ^- V ^ — ^ process with generator C^^^ defined in \2. b]) 
admits the stationary product measure: 

d S 

Poisson (A) with A := - = - . (3.12) 

7/3 
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Brownian energy process BEP(2fc) 

In this case the generator of the left reservoir is : 

. T„(24 + 4)-14, ..«^ (3,13) 

The stationary measure of this reservoir is the Gamma distribution with parameters 2k and = 2Ta. 
From the properties of the Gamma distribution one has (z) = AkT . If Ta = T\, then the process with 
generator C^^^ defined in {2.7^ admits the stationary product measure: 

Gamma (2fc,2r) with T := Ta = U . (3.14) 



4 Duality 

When the reservoirs of our boundary driven processes work at different parameters value so that 
different densities or temperatures are imposed on the two sides, the stationary measure is in general 
unknown. Remarkable exceptions are the boundary driven SEP(l), with at most one particle per site, 
for which a matrix product solution is available |DEHP| . and the case of IRW, where the product 
structure of the equilibrium invariant measure is preserved. 

An alternative approach to characterize the stationary non-equilibrium state is provided by duality. 
In section |4?T] we describe duality for the processes previously defined. Dual processes have absorbing 
boundaries at two extra sites with suitable absorbing rates depending on the parameters reservoirs. In 
general the duality functions are related to moments of the stationary distribution. In section [4.21 we 
show several applications of duality and we obtain via duality the stationary non-equilibrium measure 
of independent random walkers. 



4.1 Dual processes 

Consider the extended chain {0, 1, i, i -t- 1} obtained from the original one by adding the bounday 
sites {0, L -I- 1}. Let 77 = (771, 77^) be the configuration in the original process, we denote by 
C = (Co,'?i7 ■■■t^Lt£,l+i) £ r2'^"°' the configuration for the dual process, where the configuration space 
^Duai will be specified later. We say that {r]t)t>o and (^t)t>o are dual with duality function D(ri,£^) if 

E^[D{7jt,0] = ^i[Di7],^t)] for any t>0, (77,^ G x Qouai , (4.1) 

where denotes the expectation in the original process started from the configuration 77, whereas 
denotes the expectation in the dual process started from the configuration ^. 

Theorem 4.1. For the processes defined in Sec. 2 we have the following duality results. 



Inclusion walkers SIP(2fc). The process {r]t)t>o defined by (j2.4l) is dual to the absorbing boundaries 
process {Ct)t>o with configuration space flouai = Ng"*"^ with generator 

^'oZifiO = il-a)^i[fie-°)~fiO] (4.2) 

L-l 

+ e.(2fc+6:+i) [/(r''+') - fio] +e.+i(2fc+6) [ne'-'n no] 
1=1 

+ {13- S)^L [/(^^'^+^) - fiO] , 
with duality function 
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Exclusion walkers SEP(2j). The process (77t)t>o defined by (|2.5p is dual to the absorbing boundaries 
process {£,t)t>o with configuration space ^Duai = No x {0, 1, • • . , Sj}^ x Nq with generator 

c'oZfiO = (" + 7)ei[/(e''")-/(0] (4.4) 

L-l 

+ E - [f{e'+') - /(e)] + e.+i(2j - i^) [/(e+^-'o - /(O] 
1=1 

+ (/3 + '5)a[/(e'^''^+')-/(0] , 
wit/i duality function 

Independent random walkers IRW. The process {rit)t>o defined by (j2.6p is dual to the absorbing 
boundaries process {S.t)t>a with configuration space ^Duai = Nq^^ with generator 

^S/(e) = 76[/(e''°)-/(e)] (4.6) 

L-l 

i=l 

+ /3a [/(e^'^+') - fio] , 

wit/i duality function 

Brownian energy process BEP(2fc). The process {zt)t>Q defined by (|2.7p is ditaZ to t/ie absorbing 
boundary process {£,t)t>n with configuration space flouai ~ Nq'^^ with generator 

^DualfiO = |[/(e''°)-/(C)] (4.8) 

L-l 

+ ^e,(2fc + em)[/(r'^+') -/(?)] +e.+i(2fc+e.)[/(r+''')-/(c)] 

i=l 

+ |[/(e^'"+^)-/(e)] , 

i/ie duality function is 

D^^^iz, = (2T.)«° n zf (2T,)«-^ . (4.9) 

Theorem 14.11 can be proven by explicit computations checking that the effect of the generator of a 
process on duality functions is the same as the effect of the generator of the dual process. Hereafter 
we only include the proof of the duality property for the inclusion process, being the proofs for the 
other processes quite similar. 

Proof of Duality for the SIP(2/s). From |GKRVj we know that the free boundary inclusion process 
(i.e. the process generated by the operator Cq^^ defined in (j2.4p ) is self-dual with duality function: 
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this means that the action of Cq^^ on L'q^^(-,^) and on Dq^^{ij, •) is the same, i.e. 
thus, since Cq'^ does not act on the 0-th and L + 1-th components of ^, we have 



(4.11) 



(4.12) 



It remains to verify that the actions of the operators C^^^ and jCy/^i at the boundaries are the same 
on the duahty function. We verify this for the left boundary: 



r?i! 



a{2k + r]i) 



IV. 



771 ! 



1 + 1 -Ci)! (771-6)! 



in 
ivi - 1)! 



i5-^-(r7,0 



{2k + a 



6 



(^1 



1-6. 

SIP 



- -{0(2/0 + 771) -7(771 + 1-6)} 

- • {a(2fc + a - 1) - (7 - a)(77i + 1 - 6)} 



(??i + 1 - €1) 

(7 - a)a [i^^^^(77,ci^°) - z?^^^(77,o] = [/:g^rui^^^^(^, OKO 



(4.13) 



We have used the notations £f and C^\^\ ^ to denote the left boundary parts of the generators 
C^^^ and -Cpuai (i-e. the first line in (|2.4p . resp. (j4.2p ). By an analogous computation it is possible 
to verify that 

[Cl^^D^^^{-mn) = [C^Z^^.D^'^in.m) (4.14) 

where C^^^ and C^\!^\ ^ are the right boundary parts of the two generators. This concludes the proof 
of the duality property □ 

Remark 4.2. At this point one may wonder whether there exists a diffusion process dual to the SEP. 
All the attempts that we have done in this direction seem to suggest that this is not the case. On the 
other hand, one may extend the definition of duality at the level of the generators, i.e. we say that the 
operator C is dual to the operator Couai with duality function D{z, 77) if 



[CD{Z,.)M = [CDualD{;7Mz) 



(4.15) 



Notice that this definition does not require C and Couai to he Markov generators. Under this definition, 
it turns out that the SEP(2j) free boundary operator Cq^^ is "dual" to the differential operator defined 
in (|2.16p that has been obtained as a scaling limit of the SEP(2j). 



4.2 Moments and duality. 

In this section we provide some applications of duality. Since the dual process voids the chain, we 
show that the problem of computing stationary expectations for the original process is reduced to the 
computation of the absorption probabilities at the boundaries of the dual walkers. In partricular, we 
will see how the 7i-points correlations are related to the absorption probabilities at the extra sites 
and L + 1 of 77, dual walkers. 
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4.2.1 Stationary expectations and absorption probabilities 



In the following Proposition we provide a relation connecting the expectation of tlic duality function 
and the absorption probabilities of the dual walkers. 

Proposition 4.3. Let denote expectation with respect to the stationary measure of the processes 
defined in Section 2. Let {C{t))t>o denote the dual processes defined in Theorem \4-l\ For a given 



£, G ^Duai let Id = 'l2i=o '^'^'^ define am{£,) the absorption probabilities of the corresponding dual 
walkers initialized at ^ (i.e. S,i dual walkers start from site i), namely 

a,n{0 = mU^) - m,a+i(^) = lei - m} I {6(0) = 6, Vz = 1, . . . , L}) . (4.16) 
Then we have: in the case of the boundary driven processes SIP(2k), SEP(2j) and IRW 



m=0 



where c = ^ for SIP{2k) model, c = ^ for SEP{2j) model, c = 1 for IRW model, and where the 
densities pa and pi, are defined in Table{^ in the case of the boundary driven processes BEP{2k) 

lei 

{Diz,0)L - E (2^^-)" (2Tfc)l«l''"a™(0 ■ (4.18) 



Proof. We prove (|4.17p . Let p.L,p^,p^ be the stationary measure of the process rj with boundary 
densities pa and pt. From the definition of duality in (|4.ip and exploiting the fact that the dual 
walkers are absorbed at the boundaries, we have 

{D{v,0)l = J D{riOPL,p^,pMv) (4.19) 

t—^oc 



= hm [I?(r,,6)] 

= E (c/^-)" (CP^)'^'"™ IP? ({^o(^) = m, a+i(^) = lei - m}) , 



m=0 



where is the probability law of the dual process (^(i))t>o started at ^ at time zero, and the last 
identity follows from the formulas of the duality functions (|4.3p , (|4.5p , (|4.7p and the definitions of the 
densities given in Table 1. The proof of ()4.18p is analogous. □ 



4.2.2 Averages in the stationary state 

In this section we will see that all the boundary driven stochastic models considered so far have a 
linear density or temperature profile i.e. the expectations (rji) or (zi) with respect to the stationary 
measure is a linear function of i. This is an immediate consequence of duality since, in order to study 
the average at site i in the original process, it is enough to consider a single dual random walker 
started at i and it is an elementary fact that its absorption probabilities at the boundaries will be 
linear in i. Let us see. 

For a system of size L, the expectations (77^)^ and (zi)^ can be written, up to a factor, as the 
expectations (with respect to the stationary measures of the processes ri{t) and z{t)) of the duality 
functions D{ri,£^'^) computed in the configuration with = Sij. Furthermore, using Proposition 
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they can be explicitly found as functions of the dual absorption probabilities pi := ai{C) ^-nd 
o-oiO = ^ — Pi {o-m{0 in Proposition 14. 3p . We have 

{m)L = -{D{v,C))L^PaPr+Pb{l-Pr) i = 1, . . . , L (4.20) 

c 



for SIP, SEP and IRW, with c as in Proposition Moreover, denoting by 6a = 4:kTa and 9b = 4:kTb, 
we have 

{z,)L=2k{Diz,C))L^0a P^+9b{l-p^) i = 1, . . . , L (4.21) 

for the BEP. It remains to compute pi. 

Let Xt be the random walker moving on the chain {0, 1, . . . , L + 1} as follows. In the bulk Xt jumps 
to one of the neighboring sites with rate 1/c (with c as in Proposition (j4.3p ). whereas it is absorbed 
by the left boundary (site 0) with rate u and by the right boundary (site L + 1) with rate v. The 
values of c, u and v depend on the model, they are listed in Table ([2]) . 



System 


SIP 


SEP 


IRW 


BEP 


u 


"f — a 


7 + a 


7 


1/2 


V 


(3-5 


(3 + 5 


P 


1/2 


c 


l/2k 


l/2j 


1 


l/2fc 



Table 2: Dual processes jump rates. 

The value pi can then be interpreted as the probability for the walker Xt started at i to be absorbed 
by the left boundary, i.e. pi = P {Xoo = | Xq = i). They verify the following system of equations: 

' Po = 1 

Pi = ^P^ + ^Po 
Pi-i - 2pi + pi+i =0, i = 2, 



,L-1 



(4.22) 



PL 



-PL+l + 7^PL-l 



^ PL + l = . 

Thus Pi is a linear function of i for 1 < i < L and the solution of (|4.22p is given by: 



Pi 



L+^+^-l 

cu cv 

Hence, from ([OO)) . for SIP, SEP, and IRW we get 

Pa (L+^-z) + (*+^-l)p. 



for i = 1, . . . , L and po = 1, Pl+i = 



Ml 



(4.23) 



(4.24) 



L + ^ + ^-1 ' ' 

cu cv 

with u,v asm Table and Pa,Pb as in Table ([T]). 

Remark 4.4. Under a suitable rescaling of the constants tuning the annihilation rates at the bound- 
aries ( see Remark \2.3\) . the solutions of the exclusion and of the inclusion walkers scale to those of 
the independent walkers. 



Finally, from and (g^S]), for the BEP we get 



^,,^^?AL±1^^^±M1Z1±^ . = 1,...,L. (4.25) 
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and, by a similar computation, we find 



^'-^^ 2(L^) ^ = 1,...,L. (4.26) 



for the KMP model. 



4.2.3 Stationary product measure for the boundary driven independent walkers 

In the following proposition the stationary measure for the boundary driven IRW is obtained as an 
application of the duality property. 

Proposition 4.5. The stationary measure of the process with generator C^^^ defined in \2.6\ is the 

product measure with marginals at each site i — 1, . . . , L given by Poisson distribution with parameter 



Pa + ^ - j) + P6 (i - 1 + ^) 



(4.27) 



Proof. Since for a random variable X with Poisson distribution of parameter A the n*^ factorial 
moment is given by 'E{X{X — 1) ... (X — n + 1)) — A", to prove the proposition is enough to check 
the identity 



1=1 



To this aim consider a dual walker that starts his walk from site i S {1, . . . , L}. The probability pi of 
its ultimate absorption at site is given by 

(see (j4.23p and Table (2)). Using formula (|4.17|) and observing that the absorption probabilities of 
a total of X]i=i dual walkers, with of them initialized at site z, completely factorize because the 
walkers are independent, one has 

L 

i=l 

Inserting (j4.29p in the above formula and remembering the definition of the A^, equation (I4.28|) is 
verified and the proof of the proposition is completed. □ 



4.2.4 Duality moment functions 

It turns out from the previous section that the expectations of the duality functions D[rjt,Cj with 
respect to the probability law of the original process r^t, i.e. the "duality moment functions" 

G(?7,e,t) :=E^[Z?(77t,e)] (4.30) 

are usually some kind of moments of the original process rjt labelled by the discrete parameter ^ g 
^duai- In the case of SEP, SIP and IRW, the function G{ri,£_,t) is, up to a multiplicative constant 
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depending on ^, the ^-th factorial moment at time t when the initial value is 77. In the case of BEP, 
the function G{z,£,,t) := 'Ez[D{zt,£,)] is the standard ^-th moment. Under suitable conditions, the 
set of moments, obtained on varying the parameter ^, completely characterizes the law of the original 
process. From duality wc find that the equations for the functions G{ri,^,t) are closed and quite 
simple to write. 

Proposition 4.6. Let rjt and S,t be two dual Markov processes with duality function D{ri, ^) and let C 
and Couai be their generators, then the duality moment function G(7/,^,t) defined in (|4.30p satisfies 
the following equation: 

j^G{rj,^,t)^[CDualG{71,;t)]{0 ■ (4.31) 

Proof. For any function / = f{rj) we have 

jE^[fim)]^E,[Cf{r^t)]. (4.32) 

Given f g ^duai, applying (|4.32p to 7(77) = and using duality, namely [££>(•, ^)] (77) = 

[CDuaiD{ij, •)](0 one has 

^E^ [D{r^t,0] = [[CBu.iD{7^t, -M)] = [^Duai \D{^u •)]] (6- (4.33) 

Equation (|4.31|) follows from the definition of the function G (cfr. (|4.30|) '). □ 

Corollary 4.7. Let (■) denote expectation in the stationary state and define the "stationary duality 
moment functions" 

G{0--^{Div,0) ■ (4.34) 
Lt immediately follows from Proposition |^. 6] that G{^) satisfies the equation 

{LDuaiGm=0. (4.35) 

We will see an application of the function G in section [ 



5 Instantaneous thermalization and KMP model 

In this Section we define the boundary driven process with instantaneous thermalization. An instanta- 
neous thermalization process gives rise, for each couple of nearest neighboring sites, to an instantaneous 
redistribution of the total energy (or of the total number of particles). The class of instantaneous 
thermalization processes we consider in this paper is obtained from the non-equilibrium processes de- 
fined so far after performing a suitable "instantaneous thermalization limit" : for each bond, the total 
energy E (or the total number of particles) of that bond is redistributed according to the stationary 
measure of the original process at equilibrium on that bond, conditioned to the conservation of E. 



5.1 Thermalized models 

To start with we recall a well known instantaneous thermalization model, the KMP model (see |KMPj ). 
The KMP model is defined by considering on each bond a uniform redistribution of energy. At the 
boundaries the energy is fixed by a reservoir which imposes a Boltzmann-Gibbs exponential energy 
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distribution with different temperatures Ta and T;,. Tlie generator of the process is 

>C™/(z) - / dz[—^{f{z[,z^,...,z^)~f{z)) (5.1) 

L-1 „i 

+ X! / dx{f{zi,...,x{zi + z,+i),{l-x){zi + z^+i),...,ZL)- fiz)) 



i=l 







— ^7 — if{zi, zl-i,z'l) - f{z)) 
for any / : M. 

At the end of this section we will sec that the KMP model can be obtained as the instantaneous 
thermalization limit of the BEP(2fc) model in the particular case fc = 1/2. 

From |KMPj we know that the KMP is dual to a suitable discrete Markov process. The dual process 
C(0 = {io{t),^i{t), . . . , ^i(i), e Nq+^ describes the motion of particles in a one dimensional 
L + 2 -sites chain. The boundary sites and £,l+i are absorbing. In the bulk, for each couple of 
neighboring sites + 1) there is an instantaneous uniform redistribution of the total number of 
particles + ^i+i- The redistribution takes place whenever an exponentially distributed clock rings. 
The clocks (one for each couple (i, i + 1)) are mutually independent. The generator of this process is 
defined on functions / : Nq — > M by 

= [/(eo + 6,o,6,-.-,a+i)-/(e)] (5.2) 
+ E E [/(Co,...,e.-i,r-,6+e.+i-'^,---,a+i)-/(c)] 

i=l r=0 

+ [/(eo,...,o,a+eL+i)-/(c)] , 



and the duality function is D^^^{z,C) = n»=i T^''*' ■ 



We will see that, for each of the instantaneous thermalization processes that we are going to introduce 
there is a dual process. The dual processes are instantaneous thermalization processes themselves. 
They have absorbing boundaries and can be naturally derived by a thermalization limit from the dual 
processes of the original ones (see Section 2]). 



Thermalized Inclusion walkers Th-SIP(2fc). The instantaneous thermalization limit of the In- 
clusion process is obtained as follows. Imagine on each bond + 1) to run the SIP(2/c) dynamics 
for an infinite amount of time. Then the total number of particles on the bond will be redistributed 
according to the stationary measure on that bond, conditioned to conservation of the total number of 
particles of the bond. We consider two independent random variables rji and ry^+i distributed accord- 
ing to the stationary measure of the SIP(2A:) at the equilibrium. Thus rji and ry^+i are two Negative 
Binomial random variables of parameters 2k and p. Hence r]i + r]i-\-i is again a Negative Binomial r.v. 
with parameters 4fc and p and then the distribution of one of them, given that the sum is fixed to 
rji + = E, has a Negative Hypergeometric probability density of parameters {E, Ak — 1, 2k), i.e. 

vll^ir I E) := P(rn = v \ ry, + r;,+i = E) = ^ ^ ■ (5-3) 

I E ) 

On the other hand, the stationary distribution of the left Inclusion reservoir is the Negative Binomial 
with parameters 2k and — (resp. 4 for the right reservoir). Then the generator of the instantaneous 
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thermalization limit of the Inclusion process with reservoirs can be defined as follows 

oo 



2fc + r — l\ f a\~ f ^ — "^^ 



r J \ 1 J \ 7 

L-l iji+iji + i 

X! X! \fi^nl^■■■,Vi-l,r■,^l^ + ■n^+l~r,■qi+2,...,■nL)- fir])] {r \ i^, + -q.+i) 



It is easy to check that the thermahzed inclusion process is dual, with duality function (|4.3p . to the 
process that behaves in the bulk as the thermalized SIP(2A:) itself, and which has absorbing boundaries 
at two extra sites with absorbing rate 1. In other words the dual process is generated by: 

'Cf/,£ual/(e) = [/(eO+6,0,6,---,a+l)-/(e)] 

L-l?i+«>+i 

+ E E [/(^O' • ■ • ' ^-1' ^' + ^'^+1 - ^' ^^+2' ■ • ■ ' ^L+i) ~ /(C)] yll'^ir I + 

i=l r=0 

+ [/(Co,...,a-i,o,a + a+i)-/(0] , (5.5) 

where f = (^o, 6, • • • , Cl, 

Thermalized Exclusion walkers Th-SEP(2fc). If we take two independent random variables rji 
and 77i-|_i with Binomial distribution of parameters 2j and p, then rji + 77i-|_i is again a Binomial r.v. 
with parameters 4j and p; then the distribution of one of them, given the sum fixed to rji + rji^i = E, 
has an Hypergeometric distribution with parameters (£',4j, 2j), i.e. a probability mass function 

i^^f^ir I E) P(rn = r hi + r;^ = i?) = ^'"''l^-''^ lr<2j ■ (5.6) 

(i) 

The stationary distribution of the left Exclusion reservoir is the Binomial with parameters 2j and 7^-^ 

(resp. for the right reservoir). Then we define the generator of the instantaneous thermalization 
limit of the Exclusion process with reservoirs as follows 

CT'm = E[/(r,,„...,,,)-/(,)^ ^'-^'^ 



r=0 

L-l '/i+r/i + l 



E E \fi'ni,---,'ni-i,r.,-n, + -q,+i-r,iq,+2,----,TlL)- f{-n)\ V2^^{r\-q,+rj,+i) 



i=l r=0 



r=0 



The thermalized exclusion process is dual, with duality function (j4.5p . to the process that behaves in 
the bulk as the process itself, and which has absorbing boundaries at two extra sites with absorbing 
rate 1: 

ATDual/(e) = [/(Co + 6,0, 6,. -.,^+1)- /(e)] 

L-l C.+C. + i 

+ E E [/(a,---,6-i,^,e« + 6+i-r,c,+2,...,a)-/(6] 4r^^\^^ + ^^+^) 

z=l r=0 

+ [/(Co,...,a-i,o,a + a+i)-/(e)] , (s.s) 
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where ^ = (^o, 6, • ■ • , Ci, 



Thermalized Indepent walkers Th-IRW. Let rji and 77;+ 1 be two independent random variables 
with Poisson distribution of parameter A, then rji + rji^i is again a Poisson r.v. with parameter 2A 
then the distribution of one of them, given the sum fixed to rji + r/^+i = E, has a Binomial density of 
parameters 1/2): 

:.^«^(r I E) := F{m = r \ m + V2 = E) ^ (^^) ^ . (5.9) 

Moreover the stationary distribution of the left IRW reservoir is the Poisson with parameter ^ (resp. 
^ for the right reservoir) . Then the generator of the instantaneous thermahzation hmit of the inde- 
pendent walkers process with reservoirs is given by 

p-a/7 



Cirfiv) - E[/(^'^2,...,^L)-/(r/)] - 



,.=0 

L-l + i 

+ ^[/(^i,...,^^_i,,)-/(r;)] - — ^. (5.10) 

The thermalized independent walkers process is dual, with duality function (j4.7p . to the process that 
behaves in the bulk as the process itself, and which has absorbing boundaries at two extra sites: 

4Sai/(0 = [/(Co + 6,0, 6,..., a+i) -/(?)] 

L-l Vi+Vi+l 

+ E E [/(6,.-.,e.-i,r-,6 + e.+i-'^,6+2,...,a)-/(0] ^^'^^(ru. + c.+i) 

i=l r=0 

+ [/(eo,...,a-i,o,eL + eL+i)-/(0] , (5.11) 

where ^ = {£,0,^1, ■ ■ .,£,l,£l+i)- 

Thermalized Brownian energy process Th-BEP(2fc). We define the instantaneous thermahza- 
tion limit of the Brownian Energy process as follows. On each bond we run the BEP(2fc) for an infinite 
time. Then the energies on the bond will be redistributed according to the stationary measure on that 
bond, conditioned to the conservation of the total energy of the bond. If we take two independent 
random variables Zi and z^+i with Gamma distribution of parameters 2k and 9, then the distribution 
of one of them, given the sum fixed to Zi + Zj+i = E, has density 

J^^ zf'-^E - z^r^-^dz, 

Equivalently, the variable x ~ Zi/E has a Bcta(2A:, 2k) distribution. Denoting by V2j^^{x) the density 
of such a random variable, we can define the generator of the instantaneous thermahzation limit of 
the Brownian Energy process with reservoirs as follows 

Cf,r.f{^) = dz[ (y^)2fcr(2fc) (^i)^^'^e-^/^° {f{z[,z„ ...,zl)- f{z)) (5.13) 



p{z.\z, + z,+,^E)= ^TT— • (5-12) 



•:{zi + z^+i), (1 - x){zi + z^+i), . . . , Zi) - /(z)] V2k^{x)dx 



1 

(T,)2fer(2fc) (^^)''''^'^^^'^^ (^(^1' • ■ • ' ^^-1' - /("^)) 
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Remark 5.1. For k ~ 1/2 this reproduces the uniform redistribution rule of the KMP model on each 
bond of the bulk. The same is true for the reservoirs since the stationary distribution of the Brownian 
Energy reservoir is Gamma with parameters 2k and 9. If one takes k = 1/2, then one obtains an 
Exponential distribution with parameter 9 . 

The thermalizcd Brownian Energy process is dual, with duahty function (j4.9p to the process that 
behaves in the bulk as the thermalizcd inclusion process, and which has absorbing boundaries at two 
extra sites with absorbing rate 1. In other words the dual process is generated by: 

>cKai/(0 - [/(eo+a,o,6,---,a+i)-/(0] 

L-l?i+«, + i 

+ [/(eo,...,a-i,o,a + a+i)-/(0] , (5.14) 

where ^ = (^o, Ci, ■ ■ • , Cl, 

Remark 5.2. It is immediately seen that for k = 1/2, this gives the KMP-dual process defined in 
(EH). 

5.2 Stationary measures for L = 1. 

In this section we compute the moments of the istantaneous thermalization processes defined so far, by 
using the result obtained in Corollarv l4.7l When L = 1 there is no bulk contribution in the generator, 
since we have a site interacting with two sources. 

The L = 1 case is trivial for our original processes (SEP, SIP, IRW and BEP), because for these models 
the contributions of the baths are additive. Then the system is indeed equivalent to the system of 
one site interacting with one bath whose parameters are recombinations of the parameters of the two 
original baths. The stationary measure is, then, the stationary measure of this total bath. 
The interest of the L = 1 case for the thermalizcd processes lies in the fact that for these models 
the baths contribution arc no longer additive. Thus, even in this basic case the computation of the 
stationary measure is worth to be investigated. A result in this direction has been obtained in |BDGJL] 
(see also the remark at the end of this section). 

Thermalized Inclusion walkers Th-SIP(2fc). From the previous section we know that the ther- 
malizcd inclusion process rjt is dual to the process defined in (jS.Sp . with duality function (|4.3p . then 
the duality moment function is 

= (I>--(,.0) ^ {^f {j^f MK.) , ,5.15, 

where M(^) is the ^-th factorial moment with respect to the stationary measure of 771: 

M(er) = ij^) . (5.16) 

In order to find from Corollarv l4.7l we impose 

= /:f/,LaiG(C)=G(eo + ei,0,6)-2G(^o,ei,6) + G(6,0,a+e2) 
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This yields 



r(2fc + ei) 
2r(2fc) 



a 



"f — a 



13-5 



(5.17) 



Thermalized Exclusion walkers Th-SEP(2j). The thermalized inclusion process rjt is dual to the 
process in (|5.8p . with duality function (j4.5p . then the stationary duality moment function is 



SEP, 



r(2.7 + i-6; 



i2 



+ r(2j + l) \l3 + 6^ 

where M{S,) is the ^-th factorial moment with respect to the stationary measure of ryi: 



M(ei 



'('71-6)! 
From Corollary 14. 7[ we have 

= /:f/XaAe) = G(eo + a,o,6)-2G(eo,ei,6) + G(eo,o,6 + 6) 

r(2j + i-ei) , 



7 + a 



Co 



?2 



a 

7 + a 



2A'/(ei) ■ 



r(2j + 1) 



This yields 



r(2j + 1) 
2r(2j + i-ei) 



a 



7 + a 



(5.18) 



(5.19) 



(5.20) 



Thermalized independent random walkers Th-IRW. The thermalized IRW process rjt is dual to 
the process defined in (|5.1ip . with duality function (j4.7p . then the stationary duality moment function 

is 



M(Ci) 



(5.21) 



where A/(^i) is the ^-th factorial moment with respect to the stationary measure rji as above. To find 
M(^i) we impose 

= A'SaiG(0 = G(Co+Ci,0,6)-2G(Co,ei,6) + G(^o,0,ei+6) 

a 
7 



7/ -'''^^'^^^ \P 



This gives 



(5.22) 



Thermalized brownian energy process Th-BEP(2fc). The thermalized brownian energy process 
zt is dual to the process defined in (|5.14p . with duality function (|4.9p . then the stationary duality 
moment function is 



r(2fc) 
r(2fc + Ci) 



(2Tb)«= M(Ci) 



(5.23) 



21 



where M(^i) is now the £^1^ moment with respect to the stationary measure of zi: 

In order to find from Corollary 14.71 we impose 

- ASaiG(e)-G(Co+ei,0,e2)-2G(Co,6,6) + G(^o,0,6+6) 



This yields 



r(2fc + a 



r(2fc) 



Remark 5.3. For k 



2r(2fc) 

1/2 the above becomes 



r(2/c + Ci) 



M(6 



(2T„y^ + {2n 



(5.24) 



(5.25) 



(5.26) 



T/ie knowledge of all the moments fully describes the stationary measure of the KMP process with 1 
particle. A similar result was obtained in WDGJLf . In that paper an explicit form of the stationary 
measure for 1 particle connected to two reservoirs is given, however for a process which is slightly dif- 
ferent that the original KMP process. The difference lies at the boundaries thermalization mechanism: 
in the KMP model the first and last sites are directly thermalized by the reservoirs, in IBDGJLf the 
first and last sites share uniformly their energies with thermalized reservoirs. 



6 Correlations in the stationary state 

For some of our processes, such as for the SEP(l) model |Spo[ , the BEP(l/2) model [GKRj . the 
KMP model |BDGJLj . the covariances have been proven to be bilinear. For the boundary driven 
SEP(l), from the exact solution of the microscopic stationary state (see e.g. |DE1 IDEHPi IDLSTj ) . we 
know even more. Indeed for this process all the correlations (77^^ . . . rji^) in the stationary state are 
multilinear in the variables ii, . . . , i„, and can be explicitely computed through a recursive argument 
on n and L by a matrix method. 

The algebraic similarity of the generators for our whole class of models, that includes the SEP(l) , 
leads us to expect multilinear correlation functions. This turns out to be false in general, as we will 
see in this section. For instance bilinearity of the covariances holds only for a certain choice of the 
boundary parameters for the SEP(2j) and for the SIP(2A;) and only in the case fc = 1/4 for the BEP. 
From Proposition 14.31 multilinearity of the correlation functions is in turn implied by multilinearity 
of the absorption probabilities of the dual walkers. 

In this section we compute the two points correlations w.r. to the stationary measure, for a suitable 
choice of the parameters. We do this by direct computation, i.e. for the particle models we require 
that the generator of the process vanishes on the functions f{ri) = rjirj^. This yields a linear system 
in the variables Xi^^ := {rjirji), i < i. Analogous computations, with rji replaced by z^, are performed 
for the BEP model. 

6.1 Covariances 

The correlations in the stationary state, i.e. the expectations Xi^£ = (rjirji) with 1 < i < i < L, 
satisfy a system of L{L + l)/2 equations. The equations are quite complicated (we include them in 
the Appendix, see section ([8])) and then hard to solve directly. What we found is that, for a generic 
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choice of the boundary parameters, for none of our processes there exists a bihnear function satisfying 
them. In other words, the ansatz 

X,^e^ Aie + Bi + C£ + D for i<e and Xi^, = Ei^ + Fi + G (6.1) 

does not produce a solution for the systems in (|8.ip and (|8.ip (since the number of independent 
equations that the coefhcients in (|6.ip should satisfy is larger than 7) . But there exist some conditions 
on a,/?, 7, (5 producing an effective simplification of the systems (jS.ip . Under this conditions the 
correlations for the SIP(2A:) and for the SEP(2j') are bilinear and then explicitly computable through 
the ansatz (|6.ip . In what follows we provide such explicit bilinear forms. In order to verify their 
validity one can simply put the generic forms (j6.ip in the systems (jS.ip . find the equations that must 
be satisfied by the 7 coefficients and solve them. 

Finally, at the end of the paragraph, we will see by duality that one needs to fix fc = 1/4 in order to 
have bilinear correlations for the BEP. 

We denote by (?7i?7f)c the covariances (truncated correlations) in the stationary state of the particle 
models, namely {riirii)^ := {rjirji) — {'r]i){r]e). Replacing rji with zt one defines the covariances of the 
BEP model. 



Inclusion walkers SIP(2fc). If the parameters satisfy the condition 

7 = 2fc + a and /3 = 2fc + 5 , (6.2) 

i.e. 

cx S 

Pa = 2fc ~ a and pt, = 2k = 5 , (6-3) 

7 — a — p 

one has 

^'^^"^^^ = {L + imk{L+i) + 1) ^'"^ ' ^^-^^ 

whereas {rff) is a quadratic function of i. Notice that, under this same choice of parameters, the 
expression for the averages is simplified to 

1 

iVi) ^ Pa + [Pb - Pa)^^ ■ (6.5) 

Exclusion walkers SEP(2j). Under the choice of the parameters 

7 = 2j-a and P ^ 2j - 5 , (6.6) 

i.e. 

(1 5 
Pa = 2j — ■ — = a and pb = 2j — — = 6 , (6.7) 

a 4- 7 p + 

the two points correlations are bilinear and they are given by 



(L + l)2(2j(i + l)-l) 



{riim)c = -Trr-Tw^-yi— ^i^ — -{Pa-Pbf for i<(., (6.8) 



the variances are quadratic and the average profile becomes 



1 

{■q,) = Pa + {.Pb- Pa)j—-^- (6.9) 
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Brownian energy process BEP(l/2). In |GKRj it was studied the BEP model for k = 1/4 and it 

was found that the two points correlations are bilinear. For i < £ they are given by: 



2i{L- 
(lT3)(L + 1)2 



i^^^i)c = Zr . .\2 {Gb-eaf . (6.10) 



In this case one has the neat linear profile 



{z,) = ea + {eb-ea)-^ . (6.11) 

The result in (|6.10p can be obtained from (|6.4p and duality. Indeed, comparing (|4.2p and (|4.8p . one 
notices that the dual processes of BEP(2fc) and SIP(2fc) with 7 — a = 2fc and /3 — (5 = 2fc do coincide 
if and only if fc = 1/4. Under this choice, when the dual process is initialized from the configuration 
^ having one particle at site i and one particle at site £, equation (|4.3p becomes 

D^'^'i^hO - {2at 477,77, (2<5)«^+^ (6.12) 

and equation (|4.9p becomes 

D^^''{z,l) = (2T,)«« 4z,Z£ (2T,)«^+i . (6.13) 

Therefore, with this choice of parameters and initial conditions, the duality functions are the same if 
one identifies Ta = a = Pa and ^ S = pb and the result (jG.lOp immediately follows from (|6.4p . 



We can summarize the situation as follows. The covariances are bilinear at least in the following cases: 

(a) SEP(l), {j = 1/2) and generic a,/3,7,(5. 

(b) SEP(2j) for j G {1,3/2,2,5/2,...} and 7 + a = 2j and /3 + (S = 2j. 

(c) SIP(2fc) for fc> and 7 - a = 2k and /3 - 5 = 2k. 

(d) BEP(i), {k = 1/4) and generic Ta,Tb. 

We remark that the conditions (b),(c),(d) are those giving the neat average profile of equations (j6.5p . 
(|6.9p and (|6.1ip . i.e. those yielding exactly the densities pa and pb (resp. the temperatures Ta and 
Tb) in the proximity of the reservoirs (i.e. for i ~ and i = L + I). 

The following further properties of the covariances are observed by solving the equations in the Ap- 
pendix ([S])) on Mathematica. As the parameters are varied, the covariances are: 

(i) proportional to {pb - PaY or [Tb - TaY . 

(ii) positive for the inclusion walkers and for the Brownian energy process, negative for the exclusion 
walkers: this is related to the attractive (bosonic) interaction of the first two system, compared 
to the repulsive (fermionic) interaction of exclusion. For the proof of this property see |GRV) . 

6.2 Results for the n-points correlations 

The multilinearity of the correlations seems to be, thus, prerogative of some special cases. One may 
wonder about the multilinearity for the 3-points correlations, in the same range of parameters leading 
to bilinearity for the 2-points correlations (i.e. in the cases (b),(c),(d) above). The explicit solution of 
the n-points correlations problem becomes harder and harder as n increases and even the case n = 3 
is quite diflticult to solve exactly. 
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In this paragraph we provide the results of some numerical computations. We solved numerically the 
master equation for the invariant distribution of SEP(2j) in the cases L = 6 and j = 1/2, 1,3/2,2 
and computed the correlations (riir]j)c and {'r]i'r]jr]i)c. If ('7i?7i)c were multilinear, the differences 
di = {Tiiiji+i)c — (?7i?7i)c, i = 2,3,4,5 would be constant. The simulations seem to confirm the 
bilinearity of the covariances in the cases (a) and (b) above, and the loss of bilinearity in the other 
cases. In Fig.l (left panel) the values of di are reported for the case a = 1, 7 = 1, /3 = 1/2, 5 = 3/2: 
they are clearly constant for j = 1/2 (case (a) above) and for j ~ 1 (case (b) above) but not for 
j = 3/2 and j = 2. 

Concerning the 3-points correlations, the simulations show that the multilinearity is lost even in the 

cases where it holds for n = 2 (i.e. in the case (b)), while it is conserved for the SEP(l) with 

at most one particle per site. Figure 1 gives evidence for this phenomenon by showing that = 

iViV^Vi+i) c — {viV2Vi)c, j = 3,4, 5 are constant only for j = 1/2 (case (a) above). 

The deviation from multilinearity is in any case very small and, very likely, it is decreasing as L 

increases. 




Figure 1: Test for the multilinearity of the connected correlations {r]ir]i)c, {r]irj2'rji) c for SEP(2j), j = 
1/2 (A), i = 1 ([>), j = 3/2 (*), j = 2 (o) with a = 1, 7 = 1, ^ = 1/2, (5 = 3/2 and L = 6. In the left 
panel = ^ri\ri^+\)c - {mVi)c, « = 2, 3, 4, 5; in the right panel Ci = {iumni+i)c - {viV2ni)c, « = 3, 4, 5. 
Non constant dt or imply violation of the multilinearity. 



7 Hydrodynamic limit 

The aim of this section is to show that the large scale properties of the models studied so far can 
be obtained by a suitable adaptation of the macroscopic fluctuation theory of [BDGJLl IBDGJLH 
IBDGJL2I IBDGJL3] . In particular we verify that the macroscopic limit of the exact solutions for the 
covariances found in Section [5] does match the prediction of the macroscopic fluctuation theory (see 
jPLSi IDG) for the exclusion process with at most one particle per site) . 

7.1 Macroscopic fluctuation theory and density large deviations functional 

We briefly review the approach of the macroscopic fluctuation theory. Let us consider a one dimen- 
sional diffusive systems of linear size L in contact with two reservoirs at densities pa and pb- The 
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macroscopic fluctuation theory describes the behavior of the system in the hydrodynamic hmit in 
terms of the two quantities D{p) and cr(p), called diffusivity and mobility, defined by 



Dip) 



lim lim 

i5p^-0 L^ca op 



L {QiA+l{t)) L,p.p+Sp 



lim 

L— f oo 



(<3M+l(0)i,p.p 



where 



(7.1) 

(7.2) 
(7.3) 



In the above equation Qi^ij^i{t) is the total flow through the bond i,i + 1 in the time interval [0,i], 
while Ji,i+i(i') is the instantaneous flow at time t' . The bracket {■)L,pa.pb denotes the expectation 
with respect to the stationary state for the system of size L whose density on the left (resp. right) 
boundary is pa (resp. pb). 

From the macroscopic fluctuation theory jBDGJL2j , we know that the probability of observing a time 
dependent density and current profiles p{x, r) and j{x, r) in a macroscopic time interval [ti, T2], under 
the diffusive scaling x = i/L and r = t/L'^, is exp{~LA), where A is the action functional given 
by: 



A{{p{x,s),j{x,s)};Ti,T2) 



ds dx 



jix,s)+D{p{x,s))^£^ 



2aip{x,s)) 



Then the probability of observing a density profile p{x) in the stationary state is P{p{x)) 
where J- is the large deviation functional: 

-^(P(a^)) = , , min A{{p{x, s)J{x, s)}; ~oo,t) 

{p{x,s),j{x,s)} 



(7.4) 

g-L.F(p(x)) 

(7.5) 



with the minimum in (j7.5p taken over all the trajectories conditioned to the extreme values p{x, —00) = 
p*{x), p*{x) the typical profile, and p(x, r) = p(x). Density and current must also satisfy the continuity 
equation 

dp/dr = — dj/dx. (7.6) 



The density correlation functions in the stationary state can be obtained from the large deviation 
functional J- through the derivatives of its Legendre transform Q (see [D]): 



Q{{a{x)}) ^ sup <^ / a{x)p{x)dx - J^{{p{x)}) 

{p(x)} IJO 



Then, for large L we have 



{P{^)) 
{p{x)p{y))c 

{p{xi)p{x2) ■ . ■ p{Xk))c 



dg 



da{x) 



1 d^g 



j(2;)=0 



L da{x) da{y) 

d^g 



a(a:)=0 



1 



1 da{xi) . . . da{xk) 



a{x)=0 



(7.7) 



(7.8) 



(7.9) 
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7.2 From SEP(l) to models with constant diffusivity and quadratic mo- 
bility 

In this section we use a scaling argument to deduce the density large deviations functional of a model 
with constant diffusivity and quadratic mobility from that of the SEP(l) (cfr. also |BDGJL2] ). We 
start by recalling that for the SEP(l) one has 



D{p)^l, = 2p(l - p) , 



(7.10) 



and therefore 



■AsEP{i){{p{x, s), j{x, s)};ti,T2) = 



ds dx 



/o 4p(a;, s)(l-p(x,s)) 
then, from (|7.5p . one finds that the density large deviation functional is (see |DLSj . |BDGJL3] ) 



(7.11) 



^SEP{i){{pix)}) = sup / dx 

F(x) Jo 



p{x) loj 



F[x) 



+ {\-p{x))\oi 



1 - p{^) 

1 - F{x) 



los 



F'{x) ' 

Pa - Ph 



(7.12) 



where the supremum is taken over the monotone functions with boundary values F(0) ~ pa, F{1) — p^. 
The supremum is attained for F = Fp, monotone solution of the following differential problem: 



p{x) = F + 



F(l - F)F" 



l\2 



with 



F(0)=pa and F{l)^pb. 



(7.13) 



The connected correlation functions can be obtained by computing the derivatives of the functional 
GsEP{i) as in (j7.7p . One finds that the lowest order correlations are, for large L, 



(P(x))^^^(^) = 
(p(.)p(j/))f^^(^) = 

(p(x)p(2/)p(z))f^^(i) = 



Pail - x) + pbX 

[Pa - Pb? 



(7.14) 



L 

, {Pa - Pbf 
■ L2 



x{l - y) 
x{l-2y){l-z), 



for < a; < y < z. 



Let us now consider the generalization of (|7.10p obtained by assuming that the diffusivity is constant 
and the mobility is a quadratic function parametrized as 



D{p) = C, a{p)^2Ap{B-p), 
where A, B and C are given numbers. The action functional of this system 



A{{p{x,s),j{x,s)};Ti,T2) 



ds I dx 



AAp{x,s){B - p{x,s)) 



is related to Asep(i) through the following change of variables (cfr. jDGj ) 

C 

-^({p(2;,t), j(a;,r)};ri,T2) = —Asep{i){.{p{x, s), j{x, s)]]Cti,Ct2) 



with 



P{x,s) := —p{x,C ^s) 



and j{x,s) -.^ — j{x,C ^s) 



(7.15) 



(7.16) 

(7.17) 
(7.18) 
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The scaling (|7.17|) has been chosen among all the possible scalings connecting A and Asep(i) a-s it is 
the only one preserving the conservation law (|7.6p between p{x,s) and j{x,s). 

Then, by (|7.17p and (|7.5p it follows that the large deviation functional for the system characterized 
by (|7.15p is given by 

n{p{m = ? :FsEPii) {{B-'p{x)}) , (7.19) 



and thus 



C 



J'ilpix)}) = sup / dx 



with 



A 



p(x)\og4^ + {B-p{x))\og 



F{x) 



B - p{x) 
B - F{x) 



B\og- 



F'{x) 



Pa - Pb 



(7.20) 



Pa = Bpa, Pb = Bpb and F{x) = BF{x) , (7.21) 

where F(x) is the monotone function satisfying (j7.13p . Equivalently F is the monotone solution of 
the differential problem 



pix) ^F + 



F{B - F)F" 



Using (1721) and (fTTO)) we find 
e({a(x)}) 



with 



F(0) = Pa and F{1) = p^. 



(7.22) 



sup 



a{x)p{x) dx — J-{{p{x)}) 
sup \[ a{x)p{x)dx -^TsEP{i){{B^^p{x)})\ 

p{x)] Uo A ) 



C 



sup 

{p{x)} Uo 



1 AB 



a{x)p{x) dx - FsEP{i) {{p{x)}) 



A ' 

and, from (|7.8p and (|7.14p . we have 

dg 



^SEP(l' 



{{ABC-^a{x)}) 



(Pi^)) 
{p{x)p{y))c 
{p{x)p{y)p{z))c 



da{x) 
1 



a{x)=0 



L da{x) da{y) 

1 d'^g 



B{p{x)f^P'^^^^Pa{l-x)+p,X 

AB^ 



C 



{p{x)p{y)) 



SEP(1) 



A\ {pa - Pbf 
C 



L 



(7.23) 

(7.24) 

x{l - y) 



da{x) da{y)da{z) 



A^ 



a(x)=0 



{p{x)p{y)p{z)) 



SEP(l) 



2(^)'^^^.(l-2.)(l-z) 



and, more generally, one gets a factor i3"(j4/C)" for the n-point connected correlation function 
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7.3 Macroscopic behavior of the correlations 

With suitable choices of the parameters A, B, C we can generate the large scale limits of models that 
we have considered in the previous sections. 

Inclusion walkers SIP(2fc). For interacting particle systems, the flux across bond + 1) in a 
time interval [0,t] is given by the number of particles which jump from i to i -\- \ minus the number 
of particles which jump from i + 1 to i. i.e. 

Q^.,+i{t) - / di'[ry,+i(t')(2A: + r;,(t'))-?7«(0(2A: + ry,+i(t'))] 



= 2k / dt'[?7,+i(i')-r;.(t')] ■ (7.25) 
Jo 

As a consequence, the expectation of the flow Qi^i^i{t) in the stationary state with boundary densities 
Pa.Pb is given by 

(Q,:,»+i(i))L,P=,Pb =2k-t - {r],+i ~ ri,)L,p^,p, ■ (7.26) 
It follows, from ([Ll]) and (0211), that D{p) = 2k. 

From Section [3.2l we know that the SIP(2fc) equilibrium stationary measure at density p is the product 
of NegativeBinomial {2k, p/{p + 2k)) with a variance Var{rii) = ^^^^^p^. Using (|7.25p 

{Qh+iit))L,p,P = (2fc)2 Tdi' /'*dt"([,y,+i(i')- 77^(0] [^.+i(i")-^?.(i")])L,p,p. (7.27) 

Jo Jo 



hm I I dt' I dt"((77,+i(t')-??.(t'))('?.+i(i")-^.(i")))L,p,p 



(7.28) 



Now we have 

t-*oo t Jq jo 

= hm ^ fdt' /*dt"((,7,+i(t')-r;.(t'))(77.+i(i")-'7.(i")))L.P,P 

rOO 

= 2 di((r,,+i(0)~r7,(0))(?7,+i(0-r7.(t)))L,,,p 
Jo 

/•oo 

= 2 dt-y2{iVt+i-ni)^n[ni+iit)-rii{t)] pL,p,piv)} 
Jo ^ 

where, in the last display, PL,p,p denotes the stationary equilibrium measure. By duality, 

[n,+,{t) - ruit)] = 2fc {E, [D^'P{rit),e+')] - [D^o''^ ir,{t), O] } 

= 2k{E^.,^ [D^^^i^^m] -E?- [D^'^'ivAm} (7-29) 

where is the L-dimensional configuration (^J, . . . , ^|^) with = dij and D^^^ is the duality function 
defined in (|4.10p . Let pt{i,j) be the transition probability from the site i to the site j in the time 
interval [0, t] of a random walker on the set {1, . . . , L} moving at rate 2k, then 

E^. [D'^'P{v,m] = ^E^^- • (7.30) 

3 



29 



As a consequence (|7.28p is equal to 



2y](('7i+i - Vt)V])L^p^p ■ / dt (pt{i + l-,i)-pt{id)) 

,-1 "'0 



4 ((77i+i - ?/j)?7j+i)i„p,p • / + l,i + 1) -pt(i,i + 1)) 



OO 







OO 



= AVar{n,)- d< (pt(0,0)-ft(0,l)) (7.31) 
Jo 

where the two identities above foUow from the product character of the equihbrium measure, and from 
the fact that Ptihj) depends only on the distance \i — j\. Now the random walk pt is moving at rate 
2k, then, from the master equation we have 

2 (pt(0, 0) - ptiO, 1)) = - b*(0, 1) + ptiO, -1) - 2pt{0, 0)) = - 4 • ^ptiO, 0) . (7.32) 

2k at 

Then (|7.3ip is given by 

-2Vari7^,)- __p^(0,0)-dt = 2^.p(p + 2fc).(l-p^(0,0)) . (7.33) 

Since Poo{0, 0) vanishes as L ^> oo, we finally obtain, using (|7.2I) (t{p) = 2p{p + 2k). 
Summarizing, for the inclusion process SIP(2fc), we have 

D{p) = 2k, <j{p) = 2p{p + 2k) , (7.34) 

which implies A = —1, B = — 2fc and C = 2k. This choice produces (see (j7.24p ) the following 
correlation functions: 



(p(x)) = pa{l-x)+pb 



X 



{p{-)p{y))c = Yk '" /^^^ (i-y) 

{p{x)p{y)p{z)), = -(^l-y^i^^L-^x(l-2j/)(l-z). (7.35) 

where pa and arc the SIP(2A:) boundary densities [pa = 2fca/(7 — a) and = 2kS/{f3 — S) in 
terms of our boundary parameters). Notice that the covariances in (j7.35p do indeed agree with the 
macroscopic limit of the microscopic covariances that have been found in Section [01 fsee (|6.4p ') for a 
particular choice of the parameters. Similarly, one gets for the density large deviation functional: 



n{p{^)}) = / dx 







p{x) log ^ + {2k + p(x)) log ^, , ; + 2k log 



(7.36) 



where F ^ Fp \s the monotone solution of 



p{x) ^F+ with -F(0)=pa and F(l) = p^. (7.37) 



Exclusion walkers SEP(2j). The flux is now given by 



Qm+iW = / dt'[??.+i(t')(2i-^.(t'))-'7.:(i')(2j-r/.+i(t'))] (7.38) 







2j / dt'[rj,+^{t')-r,,{t')] . (7.39) 
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As a consequence, the expectation of Qi^i+i{t) with respect to the steady state measure reads 

{Qi,i+i{t))L.p^,p, 2j ■ t ■ {ri,+i - r]t)L,p^,p, ■ (7.40) 

Thus, from ((7J|) and we get D{p) = 2j. From Section[Owe know that the SEP(2j) stationary 

measure at density p is the product of Binomial (2j, p/2j) with a variance Varijii) = . Using 

a similar computation as for the inclusion walkers then, one can compute also the mobility, obtaining: 

D{p)=2:j, a{p)^2p{2j~p). (7.41) 

Therefore we have A =^ 1, B = 2j, C ~ 2j and, from (j7.24p . we have the following correlation 
functions: 

{p{x)) = Pa{l-x)+PbX 

{p{x)p{y)p{z)), = _(^iy^(^^i_^.^.(l_2y)(l-z) (7.42) 

where pa and pf, are the SEP(2j) boundary densities {pa — 2ja/{a + 7) and pi, = 2jS/{f5 + S) in 
terms of our boundary parameters). The second line in (|7.42p does agree with the microscopic SEP- 
covariances that have been found in Section ISTTI (see (|6.8p ) for a particular choice of the parameters. 
Moreover the density large deviation functional is given by 



T{{p{x)}) = / dx 



p{x) log + i2j - p{x)) log ^ ^ + 2j log 



Fix) '^^ " '='\2j-F{x)J ' ' "pa-Pb 

where _F = Fp is the monotone function satisfying 
Fi2j_F)F" 

in 



(7.43) 



pix)=F+ with F{Q)^pa and F{1) ^ p^. (7.44) 



Independent random walkers IRW. As observed in |DG| . the independent random walkers model, 
for which 

D{p)^l, o{p) = 2p (7.45) 

is obtained in the limit A = B^^ and C ~ 1. Under this choice, see (j7.24p . all the correlation 
functions vanish (this obviously reflects the fact that the stationary measure has a product structure, 
see Proposition 14. 5p . As S -> 00, one can see from (|7.20p that, due to the concavity of the logarithm, 
the derivative F'{x) is constant. Therefore in this limit the optimal F{x) is given by 



and one get 



dx 



p{x) lo, 



F{x) = (1 - x)pa + Xpb 

p{x) 



(1 - x)pa + Xpb 



- p{x) + (1 - x)pa + Xpb 



(7.46) 



(7.47) 



KMP model. The expectation of Qi^i+i{t) with respect to the steady state measure PL,Ta,Tt is now 
given by 



{Qi,i+iit))L,T^,n 



t- dx {[x{zi + Zi+i) - Zi] - [(1 - x){Zi + Z^+i) - 2j+l])L,T„,Tb 

Jo 

t ■ {z,+i - z,)L,T^,n (7-48) 
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then, from (|7.ip and (|4.26p . we get D{p) = 1. We know that the KMP stationary measure at 
temperature T is the product of Exponential(l/T). By a duahty argument we compute also the 
mobihty and get 

D{p) = l a{p) = 2p\ (7.49) 

The KMP model can be, then, obtained (see jDG| ) by taking the unphysical limit _B — > 0, A — ;> — 1 
with C = 1. In this limit the first three connected correlations functions (see 17.24]) are 



{p{x)) = Pa{l-x)+PbX 
[Pa - Pb? 



{p{x)p{y))c 



x{l - y) 



{p{x)p{y)p{z)), ^ -2^-P^^jJ^xil-2y){l-z), 



(7.50) 



which agree with (2.38) of [BDGJLj . Moreover the density large deviation functional that we obtain 



J-{{p{x)}) = — sup / dx 

F{x) Jo 



F{x) 



Fix) 



log 



F'jx) 

Pa - Pb 



(7.51) 



agrees with the same function computed in |BGL) . 
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8 Appendix: Equations for the two points correlations 

We provide the linear systems that must be satisfied by the two points correlation functions in the 
steady state, i.e. Xi^i = {rjirj^) with 1 < i < £ < L. In the following, equations 1),2),3) are obtained 
by letting act the generator on a couple of sites at distance larger or equal than two, equations 4), 5), 6) 
are derived from nearest-neighbouring sites, equations 7) ,8), 9) correspond to the diagonal, equation 
10) is obtained from the couple (1,L). 

Inclusion/Exclusion walkers: the equations for the inclusion walkers SIP(2fc) and for the exclusion 
walkers SEP(2j) are similar, with some relevant change of signs in the two cases; therefore we write 
them together. With the convention to use upper symbol for inclusion and lower symbol for exclusion 
in ± and =F and with the further convention that h = k ioi SIP(2fc) and h = j for SEP(2j), the 
equations read 
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1) X,^i^e + X^+i,i + X,^i,^i+X^,e+i~4Xij = for i + l<£,i>l,e<L 

2) 2h{X2,e + Xij^i+Xij+i)-{6hTa + -/)Xi^e + '2haxe = for £>2 

3) 2h{X,^L-i+X,+i^L + X,^i^L)~{6h + l3TS)X,^L + 2hSx,=0 for i<L-l 

4) hXi^^ + hX^+i^i+i + (=Fl - 4/i)Xi,,+i + + hX,^i+2 - h{xi + a;,+i) = for l<i<L~l 

5) 2/iXi4 + 2hX2a - (2(3/i ± 1) + (Ta + 7))^i.2 + 2/iXi,3 - 2hxi - 2h{l - a)x2 = 

6) 2hXL^L + 2hXL-i,L-i - (2(3/1 ± 1) + (/? T + 2/iXl_2X - 2/7.tl - 2/i(l - <5)xl-i = 

7) h{xi_i + 2xi + Xi+i) + {2h ± - AhXi^, + {2h ± = for 1 < i < L 

8) 2(2/i + (Ta + 7))^i,i - 2(2/i ± l)Xi.2 - (2/i(2a + 1) + 7 ± - 2/1x2 - 2/ia = 

9) 2(2/i + (/3 T - 2(2/i ± - {2h{25 + 1) + ± - 2/ixl-i - 2h5 = 

10) - (4/i + 7 T (5 T a + + 2/iX2,l + 2/iXi^l-i + 2/i(fai + uxl) - 

Brownian energy process BEP(2fc): the equations for the BEP(2fc) read 

1) + + + - 4Xi,f = for i + l<i, i>l,e<L 

2) 4fc(Xi,f_i + Xij+i + X2,i) - (1 + 12A:)Xi.^ + ikTa{ze) =0 for £ > 2 

3) 4fc(X,_i,i+X,+i,L+X,,i_i)-(12fc + l)X,,L+4fcTb(z,)=0 for i<L-l 

4) 2fcX,,, + 2kX,+i,,+i - 2(4fc + l)X,,,+i + 2fc(X,_i,,+i + X, ,,+2) =0 for l<i<L-l 

5) 4fc(Xi4 + ^2,2) - (12/c + 5)Xi^2 + 4fcXi,3 + 4/cTa(z2) = 

6) 4k{XL,L + Xl-i,l-i) - (12fc + 5)Xl-i,l + 4kXL-2,L + m{zL-i) = 

7) (2fc + + (2fc + -4fcXj,, = for l<i<L 

8) 2(2fc + l)Xi,2 - (4fc + l)Xi,i + 2(2fc + l)Ta{zi) = 

9) 2(2fc + 1)Xl-i,l - (4fc + 1)Xl,l + 2(2fc + l)rf,(zL) = 

10) 4fcra(zL) + ikTb{zi) - 2(1 + 4fc)Xi,L + 4fc(X2,L + = 
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